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The observation by scanning tunnelling spectroscopy (STS) of Abrikosov vortex cores in the high- 
temperature superconductor YBa 2 Cu 3 0 7 _ 5 (Y123) has revealed a robust pair of electron-hole symmetric 
states at finite subgap energy. Their interpretation remains an open question because theory predicts a 
different signature in the vortex cores, characterised by a strong zero-bias conductance peak. We present 
STS data on very homogeneous Y123 at 0.4 K revealing that the subgap features do not belong to vortices: 
they are actually observed everywhere along the surface with high spatial and energy reproducibility, even 
in the absence of magnetic field. Detailed analysis and modelling show that these states remain unpaired in 
the superconducting phase and belong to an incoherent channel which contributes to the tunnelling signal 
in parallel with the superconducting density of states. 


I. INTRODUCTION 

The interaction responsible for high-temperature super¬ 
conductivity in the cuprates has remained elusive until now. 
Various charge orders 1-5 suggest that the electron matter 
in these compounds experiences competing interactions, 
while numerous spectroscopic data support the scenario of 
preformed pairs gaining coherence at low temperature. 6-10 
The vortices offer a chance to disentangle the spectral fea¬ 
tures related to pairing from those that are unrelated. 1 In 
this context the observation of discrete states in the vortex 
cores 12-1 " has been a considerable challenge for theory, 18-25 
because the calculated tunnelling spectrum in a vortex of 
d x 2_ y 2 symmetry presents a broad continuous maximum 
centred at zero bias not observed in the experiments. 
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FIG. 1. Subgap states are observed in zero field, (a) Average 
dl/dV conductance of a 120 x 120 nm 2 area at T = 0.4 K in zero 
field. The junction resistance was 1.2 GH and the regulation volt¬ 
age 300 mV Five electron-hole symmetric spectral features and the 
background are indicated, (b) Spatial map of the superconducting 
gap A sc over the same area, (c) Histogram of the gap values 
shown in b. The colouring of bars corresponds to the color scale 
shown in b. 


Here we present high-resolution scanning tunnelling mi¬ 
croscopy (STM) measurements of YBa 2 Cu 3 0 7 _ 5 (Y123), the 
first high-temperature superconductor to reveal a pair of 
electron-hole symmetric subgap conductance peaks when 
tunnelling into a vortex core. 12 This cuprate was chosen 
to revisit the vortex core spectroscopy by STM because 
the subgap peaks are significantly stronger than analo¬ 
gous low-energy spectroscopic structures reported for other 
cuprates. 16,1 It furthermore shows spatially homogeneous 
tunnelling spectroscopy, paramount to consistently map the 
spatial extent of the subgap states in the vicinity of vortex 
cores. Our study shows that these discrete states are not 
a specific signature of the vortex cores, but are ubiquitous 
in zero field, even in spectra lacking superconducting co¬ 
herence peaks. These results reshuffle the cards for the 
vortex core description in Y123 and shed new light on its 
low-temperature non-superconducting state. 


II. RESULTS 

A. Persistence of subgap states in zero field and outside 
vortices 

We study the low-energy spectral features of Y123 in very 
homogeneous conductance maps measured along the (001) 
surface at 0.4 K. Figure 1 shows the average of 2704 spec¬ 
tra recorded over a 120 x 120 nm 2 area in zero field. The 
persistence of sharp features in this average demonstrates 
the low variability of the spectral properties across the area, 
as well as the low noise level of the instrument. On top 
of a V-shaped background with a finite zero-bias conduc¬ 
tance, which are common on as-grown Y123 surfaces, we 
identify five robust particle-hole symmetric spectroscopic 
features as labeled in Fig. 1. We shall not discuss here 
the hump, the dip, and the pseudogap (PG), but focus on 
the superconducting coherence peaks (SC) near ±17 mey 
and especially on the subgap states (SGS) at ±5 meV We 
confirm the surface homogeneity by studying the spatial 
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and statistical distributions of the superconducting gap A sc , 
defined as half the energy separation between the two SC 
peaks (Fig. lb and c). The average gap is 17.2 meV with 
a standard deviation of 1.9 meY corresponding to a ratio 
2A sc /k B T c = 4.34±0.5 where k B is the Boltzmann constant 
and T c is the critical temperature. This gap uniformity in 
Y123 contrasts with the case of Bi-based cuprates, where 
up to 50% inhomogeneity is often observed in the gap val¬ 
ues, even at optimal doping and for samples with sharp 
superconducting transitions. Most remarkable in Fig. 1 are 
the subgap states at ±5 meV Until now this structure had 
only been clearly resolved as peaks inside the vortices. 12,15 
Much weaker subgap shoulders were occasionally observed 
in zero field or outside vortices. The natural interpretation 
has been that the structure would be strongly enhanced in 
the vortex cores, suggestive of a phenomenon competing 
with superconductivity. Our measurements show that this 
is not the case: the structure is as strong in zero field as 
it is in the vortex cores. The well-defined subgap peaks in 
the average spectrum show that the energy of this structure 
does not vary appreciably along the surface. 

Applying a field of 6 T perpendicular to the Cu-0 planes, 
we observe vortices as shown in Fig. 2a. The vortices are 
imaged by mapping the ratio of the conductance at 5 and 
17 meY denoted <7 SGS /<7 SC . Outside the vortices the co¬ 
herence peaks are well developed and the ratio is typically 
0.75; this is slightly larger than the value in zero field (0.7), 
because the SC peaks are smeared due to nearby vortices. 
Inside the vortices the ratio increases up to values above 
unity as the coherence peaks get totally suppressed. The 
spatial evolution of the conductance from the vortex core to 
the inter-vortex region (Fig. 2b) shows unambiguously that 
the SGS feature is not localized inside vortices, and is the 
same structure as the one seen in zero field. The outstand¬ 
ing features of this trace are independent of direction away 
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FIG. 2. The subgap states exist inside and outside vortices, (a) 

Ratio <7 sgs /(t sc of the conductances at 5 and 17 meV measured at 
T = 0.4 K in a 6 T field over a 90 x 90 nm 2 area. The vortex cores 
appear in red and the inter-vortex regions in blue. The white arrow 
indicates a 20 nm path starting in a vortex core, (b) Conductance 
spectra along the white arrow shown in a. The spectra are coloured 
according to the conductance ratio cr s gs/ ct sc- The conductance 
scale is shown for the bottom spectrum; the others are displaced 
vertically by 0.05 nS for visibility. The distance d to the core is 
shown on the right axis. The red dashed lines mark the energies 5 
and 17 meV 


from the vortex core, although their relative amplitudes do 
depend on the proximity of neighbouring cores. 


B. Absence of correlation between subgap states and 
vortices 

Figure 2b gives the visual impression that the SGS peaks 
are reinforced in the vortex. Upon closer inspection, this 
turns out to be an illusion created by the suppression of 
the coherence peaks in the core. In fact the pseudogap, dip, 
hump, and background all remain virtually unchanged upon 
entering the vortex. The only noticeable spectral changes 
along the way are an increase of the zero-bias (ZB) conduc¬ 
tance and a suppression of the coherence peaks. The ZB 
conductance has a spatial distribution analogous to the con¬ 
ductance ratio cr SGS /cr sc displayed in Fig. 2a, although with 
a weaker contrast. This is illustrated in Fig. 3a, which shows 
a good statistical correlation between the ZB conductance 
and cr SGS /<7 sc , both increasing when approaching the vortex 
core at d = 0. As expected for a weakening superconducting 
order parameter inside vortices, we find a clear negative 
correlation between the SC peak height and cr SGS /cr sc . In 
order to quantify the amount of correlation between the 
SGS and the vortex position, we define the amplitude of the 
SGS as their height above the ZB conductance, measured 
by the difference cr SGS — <j zb . As illustrated in Fig. 3c, the 
amplitude of the subgap states is completely uncorrelated 
with the vortices. This is our central result. 


III. DISCUSSION 

The good correlation between the SC peaks and cr SGS /cr sc 
and the absence of correlation between the SGS and 
cr sGs/ cr sc imply that there is no correlation between the 
SGS and the SC peak, as we have verified in the experi¬ 
mental data. Figure 4 confirms the lack of sensitivity of 
the SGS to the superconducting coherence. A vortex-core 
spectrum (Fig. 4a) is compared with two zero-field spectra 
without coherence peaks, at a twin boundary 26 (Fig. 4b), 
and inside an extended region of the sample surface where 
no superconducting coherence was detected (Fig. 4c). All 
three spectra present the SGS feature at the same energy on 
top of different backgrounds. The PG feature appears to be 
also present in all three cases. 

STM topography reveals atomically flat terraces (rough¬ 
ness below 0.2 nm) separated by full (1.2 nm) or fractional 
(0.4-0.8 nm) unit-cell steps as well as twin boundaries. The 
absence of atomic resolution prevents a positive identifi¬ 
cation of the surface layer. But it does not preclude the 
interpretation of the tunnelling spectra in terms of intrin¬ 
sic superconducting Y123 properties as the following com¬ 
pelling evidences show. The gap defined by the coherence 
peaks satisfies the characteristic features expected for the 
superconducting gap of bulk Y123: its amplitude matches 
other bulk probes and 2A sc /k B T c is in close agreement with 
the BCS value of 4.3 for an order parameter of d x 2 _ y i sym- 
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FIG. 3. Correlation of spectral features and distance from vortex-core centre. The conductance ratio ct sgs /ct sc evaluated along the 
white arrow of Fig. 2a provides a measure of the distance d from the vortex centre. We find a strong correlation between this distance 
and (a) the zero-bias conductance and (b) the conductance at A sc , but not with (c) the amplitude of the subgap states. The colour plots 
are histograms showing the frequency of each pair of values on the same area as in Fig. 2a (white-gray: low; red-yellow; high). The 
circles and error bars indicate the average and the standard deviation of the distributions at fixed values of the conductance ratio. 


metry; the coherence peaks vanish upon warming through 
the bulk superconducting transition temperature, in the 
vicinity of twin boundaries, and when entering a vortex 
core. The very high reproducibility of the other spectral 
features as a function of surface preparation, time, applied 
magnetic field, distance from vortex cores, step edges, and 
twin boundaries makes it extremely unlikely that they re¬ 
flect electronic surface states. The latter are usually not 
as robust and one would expect modifications near step 
edges and twin boundaries. Lastly, there is a perfect match 
between the superconducting gap features measured by 
STM and point contact Andreev-Saint-James spectroscopy, 
as expected. This body of evidences unquestionably links 
A sc with Y123 bulk superconductivity and rules out elec¬ 
tronic surface states to explain the other features. We can 
also exclude that the SGS mark a proximity-induced or 
mixed d + s or d + is order parameter. In that case the SGS 
would shift to lower energy and ultimately disappear when 
approaching the vortex centres where phase coherence is 
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FIG. 4. The subgap states survive the loss of coherence. The 
thick red curves show the tunnelling conductance measured at 
T = 0.4 K in (a) a vortex core (6 T field), (b) at a twin boundary 
(zero field), and (c) in a non-superconducting region (zero field). 
The thin blue lines show a typical averaged zero-field spectrum. 
The dashed lines show that the model (1) with the red curve as 
1V 0 (5) is inconsistent with the zero-field data; a constant was added 
to match the zero-bias conductance at zero-field. All curves are 
normalized to unity at the right edge of the graph. 


lost, in clear contradiction with our data. Finally, the persis¬ 
tence of the SGS excludes their interpretation as vortex-core 
states or features of another state that would compete with 
superconductivity and pop up in the vortices. 

An obvious question is whether the non-superconducting 
spectra shown in Fig. 4 represent the normal-state DOS from 
which superconductivity emerges. Considering an isotropic 
normal state (i.e., with dispersion E k depending only on the 
modulus of k) characterized by a DOS 1V 0 (£), we find that, 
in the presence of a d -wave order parameter A cos(2ff), the 
BCS DOS can be written as 


N{E) = 


x Im 


dSNoCO — 

n 


£ + ;r + £ 


\[{E + iT) 2 - E, 2 VC E + iT) 2 - £ 2 - A 2 


( 1 ) 


The derivation is given in Supplementary Note 1. Taking 
the spectra of Fig. 4 for JV 0 (£), A = 17 me\£ and T = 
0.5 me\( we obtain the superconducting DOS shown as 
dashed lines in each panel. It is clear that the SGS are 
completely washed out, unlike in the experiment (Fig. 1). 
A feature of the normal-state DOS at energy E, indeed gets 
spread over an energy range [£, \J E, 2 + A 2 ] by the d -wave 
gap: high-energy features like the PG at E, » A remain 
sharp, but low-energy features are typically smeared over 
a range A. Illustrations can be found in Supplementary 
Figure 1. 

The picture emerging from our data is that the SGS exist 
independently of magnetic field and are not affected by 
the opening of the superconducting gap. Two possible sce¬ 
narios are compatible with the latter observation: the SGS 
either belong to an atomic layer different from the super¬ 
conducting copper oxide planes or they are confined to the 
nodal regions in momentum space in the superconducting 
layers. In the first scenario the tunnelling spectrum would 
be a superposition of two distinct contributions: one from 
non-superconducting electrons characterized by the robust 
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and homogeneous SGS as well as the ubiquitous high ZB 
conductance, another from the superconducting conden¬ 
sate responsible for the fragile coherence peaks. Indeed the 
spectrum of Fig. 1 can be perfectly reproduced by summing 
a pure d-wave spectrum and a background conductance 
similar to the spectra of Fig. 4 (see Methods). In this model 
the background forms the dominant spectral contribution 
to account for the large zero-bias and V-shaped background 
conductances systematically observed in Y123 tunnelling 
spectra. It may also explain why the vortex cores do not 
show the expected BCS spectrum (Figure 2b). Note that our 
data does show an increase of the ZB conductance in the 
cores (Figs. 3a and 4a), which may well be related to the 
predicted d-wave zero-bias anomaly that is obscured by the 
SGS. 


The second scenario assumes the SGS coexist with the 
superfluid condensate in the copper oxide planes. Then the 
SGS have to be confined to the nodal regions in momentum 
space where the superconducting gap is small. Otherwise 
they would have to shift or broaden like in Fig. 4 as the gap 
opens in the superconducting regions, in contradiction with 
the observations in Figs. 1 and 2. Supplementary Figure 2 
presents a toy model exhibiting this property. A recent study 
reports particular spectroscopic signatures from the nodal 
regions in momentum space, however without an obvious 
link to the subgap states. 28 


In either scenario it is tempting to link the SGS with the 
static charge density wave (CDW) discovered recently in 
Y123. 3,4 On one hand, the temperature scale of this phe¬ 
nomenon (~ 100 K at optimal doping, see e.g. Ref. 29) 
agrees well with the energy separation between the SGS 
peaks. On the other hand seemingly similar subgap struc¬ 
tures observed in Bi2212 vortices 14 display a modulated real- 
space pattern. 1,16,1 An explicit link between such conduc¬ 
tance patterns seen by STM and the CDW was demonstrated 
in Bi2201. 5 Preliminary data suggest similar real-space mod¬ 
ulation of the SGS in Y123, but further experiments are 
under way. 


In summary, we report the observation by STM of low- 
energy (±5 meV) particle-hole symmetric states in Y123. 
These states, seen previously only in vortex cores, are in fact 
a robust property of Y123, insensitive to the magnetic field 
and to the superconducting coherence. The subgap states 
often take the form of peaks in the tunnelling spectrum 
when the superconducting coherence peaks are absent, but 
thanks to high-resolution measurements we also see them 
as peaks even in spectra with superconducting coherence 
peaks. This new phenomenology redefines the question 
of the spectroscopy of vortices in the cuprates. Whether 
the new energy scale is a pairing scale, as its particle-hole 
symmetry suggests, whether it is connected to the Cu-0 
chains or has to do with some charge ordering, as suggested 
by similar features in Bi2212, remain open questions at this 
stage. 


IV METHODS 


Sample preparation. The experiments were conducted on 
as-grown (001) surfaces of highly pure YBa 2 Cu 3 0 7 _ 5 single 
crystals, grown in BaZr0 3 crucibles. 30 The samples were 
annealed in pure oxygen to reach an optimal doping with 
T c = 92 K. Both as-grown and chemically etched surfaces 
were investigated. As-grown surfaces were simply rinsed in 
ethanol in an ultrasonic bath before insertion into the STM. 
Chemical etching was done by immersion for 1 minute in 
a solution of 1% bromine and 99% ethanol, followed by 
rinsing in ethanol and transfer to the vacuum chamber with 
less than 5 minutes exposure to air. All data reported here 
were measured on as-grown surfaces, except the spectra 
in Figs. 4b and 4c which were measured on chemically 
etched surfaces. The spectra are independent of surface 
preparation. They are reproducible over time and after 
repeated thermal cycling and exposure to air on untreated 
surfaces, and likewise on etched surfaces if the etching is 
repeated. 


High resolution STS measurements. We used a home-built 
high vacuum variable temperature STM optimized for high 
energy resolution, with a base temperature of 0.4 K, and 
chemically etched tips made of pure iridium. Tunnelling 
spectra were acquired in zero magnetic field and in a field 
of 6 T oriented along the surface normal. The dl/dV curves, 
providing a measure of the electronic local density of states, 
were obtained by numerical differentiation of the measured 
/(V) characteristics. For measurements done on as-grown 
surfaces, the junction resistance was increased to 12 Gf2 
when moving the tip over these surfaces (as opposed to 
1.2 Gf2 in the spectroscopy mode), in order to avoid that 
contaminants stick to the tip during its motion. For the 
etched surfaces the junction resistance was kept at 1.2 Gfi 
during scanning and spectroscopy. Both procedures result 
in atomically flat surfaces allowing reproducible imaging 
and spectroscopy. 


Two-channel model for the Y123 spectrum. We assume that 
the Y123 tunnelling spectrum gets parallel contributions 
from a two-dimensional band distorted by the coupling to 
spin excitations and subject to conventional BCS d -wave 
pairing, and from a non-superconducting incoherent bath: 

dl 

— = M d JV d (eV) + M b JV b (eV). (2) 


IV d (£) is the DOS of the BCS band and M d represents 
the tunneling matrix element for this band, while M b 
and IV b (£) are the matrix element and the DOS for 
the non-superconducting bath, respectively. In the non¬ 
superconducting regions of the surface the coherent band 
turns normal with a DOS N 0 (E) and the differential conduc¬ 
tance becomes 



M d N 0 (eV) + M' b N b (eV) ss 



core 


(3) 


The matrix elements are in general different in the two types 
of regions: the constant-current STM setpoint regulation 
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FIG. 5. Two-channel model for YBa 2 Cu 3 0 7 _ 5 . (a) Best fit (solid 
blue line) of the model (4) to the tunnelling spectrum of Y123 
measured between vortices (dots), (b) Comparison of the coherent 
(green) and incoherent (red) contributions to the total spectrum. 
The numbers in parentheses give their relative contributions to the 
total spectral weight in the energy range of the figure. 


may change the relative contributions of the two channels 
when the coherent band becomes normal and its mobil¬ 
ity drops. We will assume that the spectrum measured in 
the vortex cores is a good approximation for the spectrum 
of non-superconducting regions. Eliminating N b (eV) from 
Eq. (2) by means of Eq. (3), we are left with the model 


dl 

dV 


Mb 

M' b LIdvJ core 


M'iV 0 (eV) 


+ M d N d (eV). 


(4) 


we ignore the fourth and fifth neighbour hopping ampli¬ 
tudes and adjust the chemical potential in order to reach 
a hole concentration of 0.16. The resulting parameters 
are ( t 1 ,t 2 ,t 3 ,n ) = (—281,139,-44,-356) meV Y123 is 
a bi-layer material with both bonding and anti-bonding 
bands, but there is considerable uncertainty regarding the 
precise properties of the anti-bonding band. We add an anti¬ 
bonding band via a bilayer coupling t ± and determine t± by 
least-square fitting. The two-band model for the dispersion 
is 


E = 212(005 k x + cos k y ) + 4t 2 cos k x cos k y 

t l 

+ 2t 3 (cos2k x + cos 2 fc y ) ± —(cos k x — cos fc y ) — p. (5) 

The interaction with the spin resonance is included us¬ 
ing the approach described in Ref. 32. Beside its energy 
Cl s = 41 me\( the parameters which characterize the spin 
resonance are its energy width r s and momentum width 
Aq. As the resonance is sharp in Y123 we set T s = 0. The 
resonance has a full width at half maximum 0.26 A -1 1 /a 
with a the in-plane lattice parameter of Y123. We there¬ 
fore take Aq = 1/a. The only remaining free parameter 
is a dimensionless coupling a between the Bogoliubov-de 
Gennes quasiparticles and the spin mode. We determine 
the matrix elements M b /M' b , M' d , M d , the model parame¬ 
ters tj_ and a, as well as the gap A and Dynes broadening 
T by least-square fitting of the expression (4) to the DOS 
measured in-between vortices, which is very similar to the 
zero-field spectrum shown in Fig. 1. 

The result of the fit is displayed in Fig. 5, where the fitted 
parameters are also reported. The fit yields M' d = 0, indi¬ 
cating that the coherent band is indeed invisible when it 
becomes normal; this supports an interpretation that this 
band would lie in the CuO plane below the surface and 
show up in the STM conductance only when it supercon¬ 
ducts. The model captures well the asymmetric heights 
of the coherence peaks, the relative value of the zero-bias 
and background conductances, and the dips at positive and 
negative energies. Further refinements of the model, for 
instance by allowing the adjustment of all band parameters, 
could presumably improve the agreement even further. 
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SUPPLEMENTARY NOTE 1 

One-channel models for the spectrum of Y123 in zero field. tance measured by STM on the surface of Y123 gets contribu- 
One-channel models assume that the tunnelling conduc- tions from a single electronic band of dispersion E, k = e k —/jL, 
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on top of a featureless background associated with transport 
channels of the as-grown surface and responsible for the 
large zero-bias conductance. In this frame of mind, the sub¬ 
gap peaks represent a property of the normal state which 
participates in pairing and is therefore gapped in the super¬ 
conducting state. The dispersion E, k defines a normal-state 
density of states (DOS) 
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The opening of a superconducting gap changes the DOS into 
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The quantity in brackets is the retarded single-particle 
Green’s function of a BCS superconductor with a gap func¬ 
tion A k , and T is a Dynes broadening parameter represent¬ 
ing residual impurity scattering. We restrict ourselves to 
two-dimensional models. 


Arbitrary isotropic normal state. We first consider a normal- 
state dispersion that is isotropic in the plane, E, k = c dc = 
£(k), but otherwise arbitrary. The Fermi surface is circular, 
and the details of the normal-state DOS N 0 (E) are controlled 
by the dependence of the dispersion on the modulus of 
k. Any function N 0 (E) can in principle be parametrized 
by a monotonically increasing radial dispersion q(k), by 
solving the implicit equation 2nN Q {E,)^'{k) = k{E,). Our 
goal is to express the superconducting DOS (S2) in terms 
of the arbitrary function JV 0 (E). For an s-wave gap, this is 
straightforward: 



p 00 

-i(E + ir + £) 

N S [E) = 
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— 00 
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Hence for any normal-state DOS function JV 0 (E), the one- 
channel s-wave superconducting DOS for a gap magnitude 
A can be computed by performing numerically the integral 
in (S3). If the superconducting state has d- wave symmetry, 
the situation is more complicated because of the angular 
dependence of the gap. An explicit expression can neverthe¬ 
less be obtained if the gap has the functional form A cos 2$. 
Let us write JV d (E) = — (l/7t)ImJV d (z —> E + IT) with 

C d 2 k * + g(fc) 

d __ (27t) 2 z 2 — £ 2 (k) — (Acos2i?) 2 
Working in cylindrical coordinates, we note that 
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Thus iV d (z) can be rewritten as a two-dimensional integral 
without angular momentum in the integrand, and then 
expressed in terms of the normal-state DOS: 


N d (z) 
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This leads to Eq. (1) of the main text. 

Supplementary Equation (S4) shows that if the normal- 
state DOS has sharp structures, these structures are spread 
by the opening of the d- wave gap. To see this, take 
lV 0 (O = 5{E, — E 0 ) and T = 0 + in (S4). The resulting 
superconducting DOS has three square-root divergences at 
E = E 0 and E = ±(E 2 + A 2 ) 1 ^ 2 and a square-root singu¬ 
larity at E = —E 0 . The spectral weight of the delta peak 
is spread into the intervals between these two couples of 
singularities [Supplementary Figure SI (a)]. The case of a 
normal-state DOS that is the sum of a flat background and 
two symmetric peaks at E = ±E 0 is illustrated in Supplemen¬ 
tary Figure Sl(b) for d-wave pairing and Supplementary 
Figure Sl(c) for s-wave pairing. While the peak structure 
remains visible in the s-wave case, nothing but weak fea¬ 
tures subsist in the d-wave case. In the main text it is shown 
that this one-channel model with an isotropic normal-state 
dispersion cannot explain the subgap peaks present in the 
STM spectra of Y123. A weak anisotropy of the dispersion 
will not change this conclusion qualitatively. 


Strongly anisotropic normal state with nodal peaks. The 
failure of the one-channel isotropic model to account for the 
tunnelling spectrum of Y123 suggests to devise a strongly 
anisotropic model in which the two low-energy peaks in 
the normal-state DOS would be related to features localized 
near the nodal directions in reciprocal space. These features 
would not be broadened by the opening of the d-wave gap 
and could survive in the superconducting spectrum. Such a 
model is somewhat contradictory. On the one hand, we need 
a structure with a well-defined energy to produce prominent 
peaks in the normal-state. This is usually associated with a 
localized state, a non-dispersing band, or any phenomenon 
with no dispersion in momentum space. But the absence 
of dispersion in momentum space is the quintessence of 
isotropy! 

We can escape this paradox if the two symmetric peaks in 
the normal-state DOS are due to a gap in the normal-state 
dispersion, which vanishes away from the nodal regions. 
We therefore introduce a normal-state gap Q k = E 0 6(w — 

| cos k x — cosk y \/2), where 9 is the Heaviside function. The 
parameter w controls the extension of the gapped region 
around the nodal lines, w = 0 corresponding to no gap and 
w = 1 to a fully gapped Fermi surface. The modified normal- 
state dispersion is = sign(^ fc )(^^ + 0^) 1/2 . We try this 
model on a two-dimensional square lattice with the same 
dispersion as for the two-channel model presented in the 
Methods section, Eq. (5), but without interlayer coupling 
for simplicity (tj_ = 0). The corresponding Fermi surface is 
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shown in Supplementary Figure S2(b). Using the gapped 
dispersion f fc with E 0 = 5 meV and a superconducting gap 
A k = A(cos k x — cosfc y )/2 with A = 19 mey we calculate 
the normal and superconducting DOS using Supplementary 
Equations (SI) and (S2) and r = 0.1 meV The results are 
displayed in Supplementary Figures S2(a) and S2(c). 

The model with w = 0.1 indeed has two peaks from the 
normal-state DOS surviving in the superconducting DOS. 
A new energy scale appears below E 0 . This is because the 
superconducting gap at the edges of the gapped region 




FIG. SI. Spreading of a normal-state subgap peak by the 
opening of a d-wave gap. (a) DOS calculated with Supplemen¬ 
tary Equation (S4) for a normal-state DOS jV 0 (E) = 5(E — E 0 ) 
with E 0 /A = 0.25. (b) d -wave and (c) s-wave superconducting 
DOS (blue curves) for a normal-state DOS N 0 (E ) given by a con¬ 
stant background with two symmetric peaks at E 0 = ±A/2 (red 
curves). The dashed lines correspond to a constant normal-state 
DOS JV 0 (E) = N 0 . The yellow area highlight the spectral weight 
of the peaks in the superconducting DOS, and the insets show a 
zoom of the peak region. 


is smaller than E 0 . Note that our choice of the function 
Q k with the same functional form as A fc implies that the 
edges of the gapped region are iso-contours of A k . At these 
edges we have simply A k = w A = 1.9 meV There are no 
Bogoliubov quasiparticles with energy lower than wA, such 
that the superconducting DOS is completely gapped below 
this energy. Between wA and E 0 , the momenta outside the 
gapped region contribute to the DOS, which approaches the 
w = 0 curve from below, due to the missing states from the 
gapped regions. At E 0 we have the remnant of a normal- 
state square-root singularity of weight ~ w due to the nearly 
flat dispersion of states inside the gapped region, and above 
E 0 the superconducting dispersion returns to conventional, 
and so does the DOS. 


For w = 0.5, the DOS vanishes below E 0 because w A > E 0 
and E 0 is the smallest gap on the Fermi surface. There is no 
peak at E 0 , the energy of the gapped states near the nodes 
being spread between E 0 and [E^-KwA) 2 ] 1 / 2 , where a weak 
peak remains. As seen in Supplementary Figure S2(c), the 
weight of the in-gap peak is independent of w, despite the 
fact that in the normal-state DOS this weight is proportional 
to w. 


The model of Supplementary Figure S2 is as a proof of 
principle, but we consider unlikely that it explains the sub¬ 
gap states in Y123. The latter have a spectral weight com¬ 
parable with the superconducting coherence peaks, while 
in Supplementary Figure S2 the weight of the in-gap peaks 
is tiny. More importantly, the measurements show that this 
weight remains the same in the superconducting and non¬ 
superconducting states, which is not the case in the model. 
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FIG. S2. Toy model with gapped Fermi surface in the nodal 
regions, (a) Normal (red) and superconducting (blue) DOS for 
a model with gapped normal-state dispersion along the nodal 
directions. The solid lines correspond to w = 0.1, the dashed 
ones to w = 0.5, and the dotted line shows the usual d-wave DOS 
for w = 0. (b) Fermi surface (green) and gapped region around 
the nodal lines for w = 0.1 (lightgray) and w = 0.5 (gray). The 
thickness of the green line covers all momenta contributing to the 
DOS in a. (c) Low-energy superconducting DOS for values of w 
varying between 0.05 and 0.5; the curves are offset vertically. The 
bars indicate the energies wA and [E 2 + (wA) 2 ] 1 ^ 2 . 
























